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Abstract 

Convergent expansions of the wavefunctions for the ground state and low-lying excited 
states of quantum transverse Ising systems are obtained. These expansions are employed 
to prove that the dispersion relation for a single quasi-particle has a convergent expansion. 
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1 Introduction 



The ground state and low-lying excited states of a quantum spin system, governed by a Hamil- 
tonian H, are usually studied by considering a low temperature expansion of its partition 
function e.g. [1, 4, 12, 13]. This involves an analysis of the trace of the operator, exp(— 
in the limit of large (3 (with (3 = l/ksT, ks being the Boltzmannn constant). The analysis is 
usually done by applying the method of cluster expansions to the Trotter product representa- 
tion of the operator. In contrast, Kirkwood and Thomas [5] studied the ground state of such 
systems by working directly with the Schro dinger equation, i^l^f) = -E'|\E'). They considered a 
class of Hamiltonians of quantum spin- 1/2 systems, which were small perturbations of classical 
Hamiltonians. By making a clever ansatz for the form of the ground state, they were able to 
develop a convergent expansion for the ground state which they used to show that the infinite 
volume limit of the ground state exists and has exponential decay of the truncated correlation 
functions. 

The method of Kirkwood and Thomas has been employed, in the C*-algebraic framework, to 
establish the uniqueness of translationally invariant ground states in quantum spin-1/2 systems 
[8] and has also been extended to study the ground states of models of arbitrary (though finite) 
spin, e.g. the quantum Potts model [7]. 

The approach of Kirkwood and Thomas can, however, be applied only to a restricted class 
of Hamiltonians, namely to ones that satisfy a Perron-Probenius condition. In this paper, 
we greatly extend the class of Hamiltonians to which the Kirkwood-Thomas approach can be 
applied by removing the requirement of this condition. We also simplify their method - the 
simplification being in the estimates which show that the expansions converge. The methods 
of this paper are applicable to any lattice Hamiltonian of the form 

H^-T.'^t + eY.Vx, (1) 

i X 

where Vx is a local operator that acts only on the spins at sites in X. There is a finite constant 
R such that X is only summed over subsets with g{X) < R, where g{X) is the number of bonds 
in the smallest connected set of sites containing X. Kirkwood and Thomas required that the 
matrix representing the local operator, Vx, satisfies a Perron-Probenius condition; we do not 
need such a condition. 

Throughout this paper we will work with a specific example of such a Hamiltonian, the 
Ising model in a transverse field. 

^ = -E^r+eE ^t'^h (2) 

i <ij> 

(The sum over < ij > is over all nearest neighbor sites in the lattice.) Usually perturbation 
theory in e is done with the unitarily equivalent Hamiltonian 

H^-j:<jt + eY: (3) 

i <ij> 
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However, following Kirkwood and Thomas, we use (2) in our approach. (When e = the 
groundstatc of (2), in the usual basis in which the are diagonal, is the constant wave function. 
The Kirkwood-Thomas approach is an expansion around this constant wavefunction.) An 
example of a Hamiltonian to which our methods apply, but the original approach of Kirkwood 
and Thomas does not, is the following. 

^ = -E< + ^E KM + '^'^Ml' (4) 

i <ij> 

with J 7^ 1. 

A Hamiltonian of the form 

H^-E^t + e^Vx (5) 

i X 

is unitarily equivalent to a Hamiltonian of the form (1). A more general class of Hamiltonians 
one might consider is 

H = J2Dx + eJ2V^, (6) 

X Y 

where Dx is a local diagonal operator. If the classical Hamiltonian J2x Dx has a unique ground 
state, then Tr{e~^^), and hence the ground state, may be studied by standard expansion 
methods. It would be interesting to extend the methods of Kirkwood and Thomas and this 
paper to this class of models. 

It is of interest to study, not only the properties of the ground state, but also the low-lying 
excited states of quantum spin systems. There have been a series of papers [14, 2, 10, 6, 11, 9] 
devoted to the investigation of the low-lying spectra of infinite-dimensional quantum lattice 
systems. In [14] this has been done by reducing the problem to the spectral analysis of an 
operator which is unitarily equivalent to the Hamiltonian governing the lattice system, and is 
also the generator of an infinite-dimensional Markov diffusion process. The method of cluster 
expansions has been employed to analyse the low-lying spectral properties. In particular, the 
spectrum of the Hamiltonian has been analysed in the one-particle invariant subspaces, which 
describe the elementary excitations of the ground state. In [2, 10] it has been shown that the 
one-particle excitations form a segment of absolutely continuous spectrum, which is separated 
by energy gaps from the ground state eigenvalue, and from the rest of the spectrum. 

In the second part of this paper, we study one quasi-particle excitations in a quantum spin- 
1/2 system by employing a method different from the ones mentioned above. Our method 
uses techniques, based on the Kirkwood-Thomas approach, that we develop for the ground 
state analysis. The key idea is that the lowest excited states have a definite momentum. We 
make an ansatz for the form of these states based on their momentum, and use it to develop 
a convergent expansion for them. As an application of this expansion we study the dispersion 
relation for a single quasi-particle and show that it has a convergent expansion. 

Kirkwood and Thomas also considered models with two ground states related by the global 
spin flip, e.g., the highly anisotropic XY- and Heisenberg models [5]. The simplifications and 
generalizations of their method that we introduce in this paper also apply to these types of 
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models as will be explained elsewhere [3]. We will also show that interfaces (kink states) in 
these models may by studied in much the same way that we study the excited states with one 
quasiparticle in the models of this paper. 

We consider the Hamiltonian to be restricted to finite volumes, A, which are hypercubes, 
and impose periodic boundary conditions. This restriction is essential for our treatment of 
the one quasi-particle excited states which relies on the translation symmetry of the model. 
However, for our treatment of the ground state we could consider more general A and other 
boundary conditions. 

It has been shown [11] that the spectrum of the Hamiltonian (2) above the lowest energy 
level contains a continuous part and is separated from the lowest level by an energy gap. In 
this paper we make no further analysis of the nature of the spectrum above the ground state 
eigenvalue. 

Section 2 of this paper is devoted to our version of the Kirkwood-Thomas approach to the 
ground state of quantum spin-1/2 systems, in particular to the Hamiltonian in (2). We prove 
that the wavefunction of the Hamiltonian can be expressed in the form 



where J2x hx is a function of the configuration on the lattice which can be thought of as a 
classical Hamiltonian. It is infinite range {hx 7^ for arbitrarily large X), but the size of hx 
decays exponentially as X grows. This is the form of the Kirkwood-Thomas ansatz for the 
wavefunction, but we emphasize that our approach shows that it holds even for models without 
a Perron-Frobenius condition. (Depending on the model, the function J2x hx may be complex- 
valued.) The one-particle excitations of the system are analysed in Section 3. We obtain an 
expression for the wavefunction of the lowest excited state in a subspace of fixed momentum, 
and use it to prove that the infinite volume dispersion relation for the spectral gap above the 
ground state eigenvalue has a convergent expansion. This shows that this gap persists in the 
infinite- volume limit and allows us to prove properties of the gap, e.g., that it is minimized at 
a particular value of k, depending on the sign of e. Unfortunately, the analysis of two-particle 
excitations lies outside the scope of our method. 



A classical spin configuration on the lattice is an assignment of a -|-1 or a —1 to each site in the 
lattice. So for each i e A, cTj = ±1. We will abbreviate the classical spin configuration {cri}jgA 
by (7 . For each such o" we let {(t > be the state in the Hilbert space which is the tensor product 
of a spin up state at each site with o"j = +1 and a spin down state at each site with ai = — 1. 
Thus |cr > is an eigenstate of all the with a^\a >= ai\a >. Any state > can be written 
in terms of this basis: 



^jj{a) = exp(^/ix((7)) 



(7) 



X 



2 The ground state 




(8) 



a 
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where i^ia) is a complex- valued function on the spin configurations a. 

The form of the Hamiltonian in equation (3) seems natural for perturbation theory in e 
since the Hamiltonian is diagonal when e = 0. However, following Kirkwood and Thomas, we 
study a unitarily equivalent Hamiltonian 

^--E^r + ^E ^M- (9) 

i <ij> 

For a single site the ground state of —cr^ is the vector | + 1 > +| — 1 >. Thus the unnormalized 
ground state of the Hamiltonian (9) with e = is given by (8) with ip{a) = 1 for all a. 
For a subset X of the lattice A, let 

a{X) = n ^i, 

and wc use the convention that (t(0) = 1. Let /(cr) be a complex- valued function of the cr's. 
Then /(cr) can be written in a unique way as 

f{a)^^g{X)a{X), 

X 

where the g{X) are complex numbers and X is summed over all subsets of A (including A itself 
and the empty set). The Fourier coefficients g{X) are given by 

g{X)^2-\^\^a{X)f{a), 

a 

where |A| is the number of sites in A and the sum is over all 2'^' spin configurations on A. [In 
general, let \X\ denote the number of sites in any subset X of the lattice.] 

The heart of the Kirkwood Thomas method is the following. We expand the ground state 
with respect to the basis as in eq. (8) and assume that ip{a) can be written in the form 

V^(<7)=exp[-V^(X)<7(X)]. (10) 

^ X 

(The factor of — | is introduced for later convenience.) Kirkwood and Thomas use a Perron- 
Frobenius argument to show that ip{a) > for all a, and so the ground state can be written 
in this form with real g{X). The Perron- Probenius argument requires that the off-diagonal 
terms in the Hamiltonian are all non-positive, and so it appears that the Kirkwood Thomas 
approach only applies to a very restricted class of models. This is not the case. We will show 
how to eliminate the Perron-Frobenius argument entirely. We proceed in two steps. First, we 
will make the ansatz (10) and prove that there is an eigenstate of this form. Later we will give 
a short argument to show that this eigenstate must in fact be the ground state. 
Consider the eigenvalue equation 

i/|^>=Eo|^>. (11) 



5 



The operator erf (j| is diagonal, so 

^^^|E^(^)I^ >=Y.^^^3^i^)\^ > ■ (12) 
a a 

The operator uf just flips the spin at site i, i.e., a^\a >— |(7^*) >, where crW is the spin 
configuration cr but with (jj replaced by — (Tj. Hence 

>= EV'Wk^^^ >=J:^P{a^'^)\a > . (13) 

<T (7 (T 

The last equality follows by a change of variables in the sum. 

Thus if we use (8) in the Schrodinger equation (11) and pick out the coefficient of \a >, 
then for each spin configuration a we have 

_^^(^«) + , ^ a,aj^l;{a) ^ Eo^Pia). (14) 

i <ij> 

Dividing both sides by ip{a), we have 

Now a^^{X) is just a(X) ii i ^ X and -a(X) iiieX. Thus 

V^(a«) = exp[-i g{X)a{X) + l ^ ^(X)a(X)], 

and so we have 

^^ = eMj: 9{X)a{X)]. 
Thus the Schrodinger equation is now 

-Eexp[ E 9iX)a{X)] + eY: '^i'^j = ^o- (16) 

« X:ieX <ij> 

This equation must hold for every configuration cr. This is the magic step in the Kirkwood 
Thomas method because of the cancellation between terms in ■'tp(a) and ip{a^^^). As |A| oo, 
the sum J2x 9{X)'^{X) will diverge, but the sum J2x:iex 9{X)'^{X) th^-t appears in the above 
equation will remain bounded. We will refer to equation (16) as the Kirkwood Thomas equation. 
We define 

exp(^)(x) = - 1 -X. (17) 

and rewrite (16) as 



-E^xp^'M E 9{X)a{X))+eY: cr.a, = Eo + |A| + E E 9{X)a{X) 

i X:ieX 

Eo + \A\+Y,\X\g{X)a{X) (18) 



So the Kirkwood Thomas equation becomes 

^\X\g{X)a{X) + Eo + \M^-^e^p^'^( g{X)a{X)) + e Y: '^i'^r (19) 

X i ^ X:i£X ^ <ij> 

Expanding exp^^^ in a power series we have 

Y\X\g{X)a{X) + Eo + \K\ = -EEri Y 9iX)a{X)) + 6 ^ 

X i n=2 ^X:ieX ' <i3> 

oo -| n 

= - E E 3 E n 9{X,)a{X,) + 6 E 

i n=2 "■■ Xi,X2,-,X„:ieXj j=l <ij> 

(20) 

The notation : i E Xj means that i belongs to Xj for all j = l,2,---,n. Since = 1, 
a{X)a{Y) = a{XAY) where the symmetric difference XAY of X and Y is defined by XAY = 
(X U r) \ (X n y). Thus n"=i cr{Xj) = (7(Xi A • • • AX„). Eq. (20) implies that the coefficients 
of cr{X) on both sides of the equation must be equal. So for X 7^ we have 

^(X) = -^EE^ E 9iXi)g{X2)---g{Xr.) + ^^elnn{X), (21) 

where lnn(-^) is 1 if X consists of two nearest neighbor sites and is otherwise. 
Similarly, for X = 0, we obtain the equation 

^o + |A| = -EE A ^ g{X,)g{X,)...g{X^), (22) 

i n=2 Xi,X2,--,Xn:ieXj, 
X]^A---AX„=0 

since (t(0) = 1. 

We let g denote the collection of coefficients {g{X) : X G A, X ^ 0}, and think of eq. (21) 
as a fixed point equation, g = F{g). We define a norm by 

ll^ll = E 1^(^)11^1(^1^)""'^''^ (23) 

X3i 

where w{X) is the number of bonds in the smallest connected set of bonds that contains X. 
(We make the convention that a bond contains its endpoints.) If X just contains the single site 
i, then w{X) = 0. With periodic boundary conditions we need only consider g{X) which are 
translation invariant, and so the above quantity does not depend on i. With other boundary 
conditions we would need to define the norm to be the supremum over i of the above. The 
constant M will be chosen later. It will only depend on the lattice. We will prove that the 
fixed point equation has a solution in the Banach space defined by this norm if |e| is sufficiently 
small by applying the contraction mapping theorem. 
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Theorem 1 There exists a constant M > which depends only on the number of dimensions 
such that if\e\M < 1, then the fixed point equation (21) has a solution g, and \ \g\\ < c for some 
constant c which depends only on the lattice. 

Proof: We will prove that F is a contraction on a small ball about the origin, and that it maps 
this ball back into itself. The contraction mapping theorem will then imply that F has a fixed 
point in this ball. We denote the radius of the ball by 5 and let 



M 



(24) 



(The factor of Ad is for the lattice Z'^. In general this factor is twice the coordination number 
of the lattice.) For the sake of concreteness, we will prove F is a contraction with constant 1/2, 
but there is nothing special about the choice of 1/2. So we need to show 



\\F{g)-F{g')\\<]^\\g-g'\\ for 1 1^| |, 1 | < 5, 



and 



||F(^)||<5 for \\g\\<5. 
The proof of (25) proceeds as follows. 

OO -j^ 



(25) 
(26) 



n9)-n9')\\ < E E rr E • ■■9{Xn)-g'{x,) . . •^'(X„)|(|6|M)-(^^^-^^"). 



i n=2 Xj^,-,x„:ieXj, 



(27) 

If e A . . . AX„, then is in at least one Xj. Using the symmetry under permutations 
of the Xj, we can take e Xi at the cost of a factor of n. We claim that if i e Xj for 
j — 1,2, ■ ■ ■ ,n, then 

n 

(28) 



w(XiA---AXn)<J2^(^j)- 

i=i 



To prove the claim, let Cj be connected sets of bonds such that Xj C Cj and \Cj\ = w{Xj). 
Define C = U^^^Cj. Since all the Xj contain i, all the Cj contain i. So C is connected. Clearly, 
XiA ■ ■ ■ AXn c C. So 



w{XiA ■ ■ ■ AXn) <\C\<J2 = E^W 

j=i j=i 



(29) 



which proves the claim (28). So 

m9)-F{9')\\ 



<EE 



n=2 



[n-iy. 



E 



Xi,---,X„:i^Xj,0&Xi 



n gi^j) - n 9'{x,) 
j=i j=i 



{\e\M) 



-w{Xl)-HX2y—w{Xn) 
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oo -j^ n / fc— 1 



<E7r3Tv E E E E[EmM9ix,)-9'ix,)\i n I^WI) 

n=2 V"- Xi30 ieXi X2...X„3i fe=l ^ j=l ' ^j=fe+l ^ 

oo -| / \ 

<E7-3wE E I^W-^WKI^I^)""'^''^^ E ni5'WI(l^l^)""^''^M 



n=2 V" Xi90 iGXi X2...X„9i ^j=2 

oo -1 / k—1 

+ E7-3W E E I^WKI^I^)""'^''^^ E E (ni^W 1(1^1^)""'^''^^ 



X |^(X,)-^'(X,)|(|6|M)-(^'=) n |^'(^,)|(|6|M)-(^^) 



oo 1 

<II^-^1IE 



i^'ir'+Eii^ii'"'ii^"'"-' 



k=2 



(30) 



where we have bounded various sums by \ \g — g'\\ ,\\g\\ and in obvious ways. Since 
and \ \g'\ \ are both bounded by 5, it follows from (30) that 



\F{g)-F{g')\\<K\\g-g'\\, (31) 



where 



n 



^ = E7-^rTT^"-' = e^-l + <^e^ (32) 

n=2 \"' 

Recalling that 6 = 4d/M, if M is large enough then K <l/2. 

To prove (26), we use (25) with g' = 0. Using (21) to compute F{0), it follows that 

\\F{g)\\ < ||i^(^?)-i^(0)|| + ||F(0)||<^||^||+Ek|lnn(X)(|6|M)-W 

1 S 
< -\\g\\ + 2d\e\{\e\M)-^ < - + 2dM-^ ^ 6. (33) 

■ 

At this point all we have shown is that we can find an eigenfunction of the form (10). 
Kirkwood and Thomas show that this eigenfunction is the ground state by a Perron-Probenius 
argument. This limits their method to models in which the off-diagonal matrix elements of 
the Hamiltonian are non-positive. We give a different argument which does not impose any 
restrictions on the signs of the matrix elements of the Hamiltonian. 

When e = 0, we have g = 0, and we know that our eigenfunction is the ground state. Since 
we are working on a finite lattice, we are considering a finite dimensional eigenvalue problem. 
So the only way our eigenvalue can cease being the ground state energy is if it crosses another 
eigenvalue. So if we can show that our eigenfunction is nondegenerate for all e with |e| < 1/M, 
then it follows that our eigenfunction is the ground state for all such e. Suppose that there 
is some e such that the degeneracy of the eigenvalue of our eigenfunction is greater than one. 
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Recall that our eigenfunction is denoted by ip{a). Let 0((t) denote a different eigenfunction 
with the same eigenvalue. Define 

ijjM^tjj{a) + a(t)ia). (34) 

So ipa is also an eigenfunction with the same eigenvalue. Since V'(cr) > for all a, and there 
are finitely many a, if a is sufficiently small, then ipa{<^) > for all a. Thus it can be written 
in the form 

M'^)^eM-\j29a{X)a{X)]. (35) 
^ X 

Moreover, as a ^ 0, US' — fifall 0. Since ip^ is an eigenfunction, satisfies the fixed point 
equation. For sufficiently small a, is in the ball of radius 5 about the origin. But this 
contradicts the uniqueness of the solution of the fixed point equation in this ball. Thus the 
eigenvalue is nondegenerate, and so the eigenfunction we have constructed is indeed the ground 
state. In the preceding argument, how small a must be depends on A. This is not a problem. 
For each A we need only find one value of a which leads to a contradiction and so implies that 
the eigenfunction is not degenerate. This argument works for all e with |e| < 1/M, and M does 
not depend on A. 

Finally, we consider the question of convergence to the infinite volume limit. One approach 
would be to show that the coefficients g{X) have convergent expansions in e and that these 
expansions agree up to order where L is the length of the smallest side in A. We will sketch 
a different approach based on the fixed point equation. It takes advantage of the fact that if 
we have a g which is near the origin and very close to being a fixed point, then it must be very 
close to the true fixed point. 

We start by asking where the vohime dependence is in the fixed point equation (21). The 
Xi in this equation must be subsets of the volume A, and the notion of nearest neighbor in 
the term lnn(-^) in (21) depends on the volume because of the periodic boundary conditions. 
To make the volume dependence of F[g) explicit, we write it as Fj^{g). We can define a fixed 
point equation which should give the infinite volume limit directly. To do this, in (21) we let 
the Xi be subsets of Z'^ and interpret lnn(-^) to mean nearest neighbors in Z"^ in the usual 
sense. We denote the resulting function by F^{g). Our proof of Theorem 1 also shows that 
there is a small neighborhood of the origin in which the fixed point equation F^[g) = g has a 
unique solution. 

Let gA and goo denote the fixed points of Fa and F^q. We want to show that as the volume 
A converges to Z*^, g\ converges to g^o- Note that g\ and g^o are defined for different collections 
of subsets, so it is not immediately clear in what sense g^ should converges to g^o- In fact, if 
a set X is near the boundary, g\{X) and goQ{X) need not be close. Consider one dimension 
and let A = {1, 2, • ■ ■ , L}, and let X = {l,L}. Because of the periodic boundary conditions, 
the two sites in X are actually nearest neighbors. So gA{X) is of order e. But goo{X) will be 
of order e^. 

To compare g\ and g^o correctly, we proceed as follows. We take our volume to be centered 
about the origin and let L be the length of the shortest side. So the closest sides are a distance 
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L/2 from the origin. Recall that w{X) is the number of bonds in the smallest set of bonds 
which contains X. We claim that 

E \9k{X) - g^{X)\ (36) 

X30:w(X)<L/4 

converges to zero exponentially fast as A converges to Z'^. Note that the terms we are neglecting, 

E I^^aWI (37) 

XcA:0eX,u;(X)>L/4 

and 

E \9oo{X)l (38) 

XcZ'i-.QeX ,w(X)>L / i 

are both of order (|e|M)^/^ by our bounds on the solution of the fixed point equations. 

Our strategy is to show that (?a is almost a fixed point of Foo- This then implies that giA is 
close to the exact fixed point 5^00 of Fqq. To show that gx is almost a fixed point of Foo we need 
to estimate 

i^oo(^A) - = i^oo(^A) - i^A(^A). (39) 

However, g\ is only defined on subsets of A, so we cannot plug directly into F^o- We define a 
function (7^ which is defined for all X C .Z^*^ as follows. For sets X with e X and w{X) < L/A, 
we let g'x{X) = gx{X). We then extend the definition to sets which are translates of such X 
by taking g'jt^{X + t) = g'j^{X). (These translations are done in Z'^.) Finally we define g'\{X) to 
be zero for the sets for which it is not yet defined. Note that there is a consistency issue here. 
If we have two sets X and Y which both contain 0, have w{X),w{Y) < L/4 and are translates 
of each other in Z'^, then we need to be sure gA{X) — 5'a(5^)- These conditions on X and Y 
ensure that they are also translates in A, so this follows from the translation invariance of g\. 

We will now show that gf^ is close to g^ by showing that g'j^^ is almost a fixed point of F^. 
We will estimate 

Eii^oo(^;)(^)-^;(^)i (40) 

where Foo{g'/<J{X) is the coefficient of (t{X) in Foo{g'j^). Note that we are now using a norm 
without the weighting factor of {\€\M)~^^^\ We divide the terms in this sum into two classes, 
those with w{X) < L/4 and those with w{X) > L/4. For the latter, g'/^{X) = and 

E l^oo(^;)(^)| < (|6|M)^/^ E |Foo(^;)(^)|(|6|M)-W = {\e\M)^/'\\F^{g'^)\\ 

X30:w{X)>L/4 X30:w(X)>L/4 

(41) 

where ||-Foo(5'a)II the original norm with the exponential weight. Now consider an X con- 
taining for which w{X) < L/4. Then g'/^{X) = gA{X) = Fa(5'a)(-^)- So we need to compare 
FoQ{g'/^){X) and -Fa(5'a)(A"). Consider a term Xi, ■ ■ ■ , X^ in F^o with J2i'^{Xi) < L/4 and 
XiA ■ ■ ■ AXn = X, with X containing the origin. If w(Xj) < L/4, then X^ has a translate 
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Xi-\-t with e Xj + i and w{Xi-\-t) < L/4. So we have g'{^{Xi) — g^iXi), and the contribution 
of this term to Foo{gA){X) will be the same as the contribution of this term to F\[g\)[X). 
So the difference between Fooig'/JiX) and F\[g\){X) comes from terms with J2iw{Xi) > L/4. 
Their total contribution to (40) is of order (|e|M)^/^. Combining this with (41), we have shown 
(40) is of order (|e|M)^/^. If we apply the contraction mapping theorem with the norm without 
the weighting factor (|e|M)"'(^\ this imphes 

- 9oo{X)\ (42) 
is of order (\e\M)^^^. And so (36) converges to zero like (|e|M)^/^ in the infinite volume limit. 

3 Excited states 

In this section we study the low-lying excited states with one quasiparticle for the model defined 
by (9) when e is small. When e is 0, these states for the Hamiltonian (3) have all spins pointing 
in the ^-direction except for one spin pointing in the opposite direction. Thus the energy of 
the first excited state is |A|-fold degenerate (|A| is the number of sites). When e is small but 
not zero, this degenerate eigenvalue should open up into a band of continuous spectrum in 
the infinite volume limit. The goal is to get a convergent expansion for the eigenfunctions 
corresponding to this continuous spectrum. This looks impossible since there is no gap around 
these eigenvalues. However, these excited states are uniquely determined by their momentum. 
If Ti denotes translation by /, then an eigenstate with momentum k satisfies 

Ti\^>^ e-^'^'l* > . (43) 

(For d > 1, k and I are vectors, and kl will denote their dot product.) If you look in the 
subspace with momentum k, then the first excited state has an isolated eigenvalue. We will 
show that by building into the wave function the constraint that it has momentum k, one can 
get a convergent expansion for the eigenfunction. One of the consequences of our expansion is 
the following theorem about the energy dispersion relation. 
We will always work on a hyper-rectangle: 

A = {(zi, ■■■,id):l<ij< LjJ - 1, 2, • • • d} (44) 

where is the length in the ith direction. The possible momenta are given by k — {ki, ■ • ■ ,kd) 
where ki = and = 0, 1, Lj — 1. We denote the ground state energy for a volume A 
by Eq. We let Ei{k) denote the lowest eigenvalue in the subspace of momentum A: if A; 7^ and 
the next to lowest eigenvalue if A; = 0. (We think of this as the energy of the state with one 
quasiparticle with momentum k.) These eigenvalues depend on the volume A. When we wish 
to emphasize this dependence, we will write them as Eq and E^{k). The dispersion relation is 

^\k) = E^{k) - E^ (45) 
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When e = 0, A^{k) — 2 for all k. For small nonzero e we expect A^(A;) to equal 2 plus a 
function that is of order e. The following theorem essentially says that this function has an 
infinite volume limit, and this limit has a convergent expansion. We state the theorem in terms 
of the Fourier coefficients of A^(A;) rather than the function itself since A^{k) is defined for 
different values of k for different A. 

Theorem 2 There exists an eg > such that for all \e\ < eg the following is true. Write A^{k) 
as a Fourier series: 

A^{k)^2 + J2eje"''. (46) 

There are coefficients e^, defined for all s G Z"^, such that for any sequence of finite volumes 
which are hyper-rectangles and converge to Z'^ , the coefficients e^ converge uniformly to e^. 
The coefficients Cq and cq are at least first order in e, and e^ and Cg are at least of order e'*' in 
the sense that there is a constant c such that 

|e,| < (c|e|)l< (47) 

Thus the infinite volume dispersion relation may be defined by 

A{k) =2+J2 e^e''^ (48) 



Remark: The bounds on the eg are only upper bounds. Some of the may be higher order 
in e than the theorem indicates. 

We start the proof by considering what states with momentum k look like. A function ip{a) 
has momentum k if 

Tjilj{a)^e-'''^tP{a). (49) 
If X is any subset of the lattice, then we can form such a function by letting 

</'x,fc((T) = 5:e*'=V(X + 0, (50) 
I 

where X + 1 is the set of sites {i + l : i E X}. The above sum is over sites in A. All subsequent 
sums in this section will be over sites in A unless otherwise indicated. We claim that these 
functions span the subspace of momentum k. It suffices to show that if we consider all k, then 
these functions span the entire state space. Since 

^0x,fe(<7) = |A|(j(X), (51) 

k 

the span of the functions (f)x,k{<^) includes the functions cr(X) for all X. The latter form a basis 
for the state space, so the (j)x,k span the state space. 
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While these functions span the subspace of momentum A;, they are not hnearly independent. 
For some choices of X and k, will be zero. For example, if X = or A, then the function 

is zero unless A; = 0. Furthermore, 

(l>x+t,M^e-'''(f)x,M- (52) 

So if two X's are translates of each other, then the corresponding functions will be the same 
up to a multiplicative constant. Define X and Y to be translationally equivalent ii X = Y + n 
for some n E A. Then the subsets of A can be partitioned into equivalence classes. Pick one 
set from each equivalence class and let X be the resulting collection of sets. Every subset of A 
can be written in the form X + n for some X E X. While X is uniquely determined, n is not. 
(For example, consider X = or A.) If we consider the collection 0x,ik for X E X, then this 
collection spans the subspace of momentum k. 

Now suppose that for each k we have an eigenstate ipkic) with momentum k. Let 
denote the groundstate wavefunction, (10), obtained in Section 2. Then '^k[a)/'^{a) also has 
momentum k. So it can be written in the form 

^= Y^c{X,k)M<y) (53) 

for some coefficients c(X, /c), which depend on k. For each X, c(X, k) has a Fourier series 

c(X, k)^Y. e"'*"" e(X, n) . (54) 

n 

(The sum over n is over the sites in A.) Then 

V'.(a) = V(^) E Ee"''"e(X,n)(/)x,.(^7). (55) 

XeX n 



Using (52) this becomes 



^k{(y) = ^(o-) E E ^(^' ^) (i^x+nM = ^(o-) E ^(^) '/'^,fc(^)' (56) 

X€X n X 

where the coefficients e{X) are given by 

e{X)= E e(y,n). (57) 

YeX,n-Y+n=X 

We now have 

V'fe(<7) = V'(<7) Ee^'^ Ee(^) '^(^ + 0- (58) 

I X 

If we define 

0(a) = Ee(^)'^(^) (59) 
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then 



(60) 



Note that the only momentum dependence in the above is in the exphcit factor of e^^K 

When e = 0, the lowest excited states with momentum k of the Hamiltonian H defined in 
(9) are 

-Y.^'^' ^'i- (61) 



In the above formulation this corresponds to c({0}, k) = \ for all k and c(X, A;) = for all k 
and all X e A' other than {0}. This is equivalent to e({0}) = 1 and e(X) = for all X other 
than {0}. For nonzero e we can normalize the eigenstates so that c({0}, k) — 1 for all k. This 
is equivalent to taking e({0}) = 1 and e({/}) = for / 7^ 0. We will show that the other e(X) 
arc small and localized near when e is small. As we showed in the last section, the ground 
state wave function V'((t) can be written as 

^{a)^e^v[-\Y.9{X)a{X)] 
^ X 

The ground state is translationally invariant, so g{X + /) = g{X). 
The action of the Hamiltonian on the wavefunction (58) is 



mk){cj) 



mi: 



Akl 



-Eexp(E^(mn)[ E e{X)a{X + l)- e{X)a{X + l) 

j YBj X^j-l X3j-l 

+e ^ a,a,Ie{X)a{X + l) 

<ij> X 



The above needs to equal 

Ei{k)M(^) = E,{k)ilj{a)Ye"''Ye{X)a{X + 0- 

I X 

Cancelling the common factor of ^(cr), the Schrodinger equation for the excited state becomes 



E 



leME 9iY)a{Y))[ e(X)a(X + l)- ^ e(X)a(X + l)] 

j Y5j X^j-l X3j-l 



+e Y a,a,Y<X)a{X + 1) - Ei{k)Y<X)a{X + 1) 

<ij> X X 



= 0. 



Recall that the ground state wave function satisfies 



-Eexp(E^(^M^)) + ^E (^i^^J-Eo. 

j Y3j <ij> 



(62) 



(63) 
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Multiplying this equation by J2i ^^^^ Y.x e{X)a{X + 1) and subtracting the result from eq. (62), 
we have 

2Y.e''Y.^MT.9{Y)a{Y)) e{X)a{X + l) 

I j YBj X3j-l 

= Y: e'''[E^{k) - Eo] E e{X)a{X + I). (64) 

I X 

Define h(Y) by 

eMj:9iY)a{Y)) ^l + Yh{Y)a{Y). (65) 

Y90 Y 

Explicitly, 

OO 1 

n=l Y-i, --,Yr^:QeYj, 

Note that h{Y) is localized near the origin. Using the translation invariance of the g{Y), we 
have 

exp(E 9{Y)a{Y)) = 1 + E h{Y)a{Y + j). (67) 

Y3j Y 

Inserting this in (64) 

2Ee^'T[l + EMn^(^+j)] E e{X)a{X + l) 

I j Y X3j-l 

= Ye'^'lEiik) - Eo]Y<X)a{X + 1). (68) 
; X 

The change of variables j ^ j + 1 changes this to 

2 E e^'' E[l + E h{Y)a{Y + j + /)] E e(X)a(X + /) 

I j Y X3j 

= E e'''[E^{k) - Eo] E e{X)a{X + I). (69) 

I X 

Using 

E E e{X)a{X + = E l^|e(X)<7(X + /) (70) 

j XBj X 

this becomes 

2 E 6^^=' E \X\e{X)aiX + + 2 E e''' E E E Mne(^)<7([(y + i) AX] + 
i X I Y X jex 

= Ee'''[^iW-^o]E<^M^ + 0- (71) 
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We can write Ei{k) — Eq as a Fourier series in k. Since this quantity equals 2 when e = 0, 
it is natural to write it in the form 

E,{k)-Eo = 2 + J2e-"'"'em. (72) 

Inserting this in (71) and making the change of variables Z — > Z + m and then X ^ X — m in 
the second term, the right side of (71) becomes 

2 J2 e*'' E e(^)a(X + + e^'^e^ E <^ ' + (73) 

IX l,m X 

Thus eq. (71) is of the form 

Ee^'T/(^M^ + = O (74) 

I X 

where 

f{X) = 2\X\e{X)+2 Yl h{Y)e{Z)-2e{X)-Y,eme{X-m). (75) 

Y,Z,jeZ:{Y+j)AZ=X rn 

Obviously, if f{X) ~ for all X, then eq. (74) holds. We will argue that the converse is true. 
We rewrite (74) as 

Ee^'T/(^-0^(^)=O. (76) 

I X 

This must hold for all ex, so for all X and k we must have 

Y.e'^'f{X-l)^Q. (77) 
I 

Fix an X and regard f{X — /) as a function of /. Since the above equation holds for all A;, we 
must have /(X — /) = for all /. This proves the converse, and so eq. (74) is equivalent to 
/(X) = for all X. 

We now assume that e({0}) — 1 and e({/}) = for alH 7^ 0. A priori there is no reason that 
a solution with these properties must exist, but we will show that it does. As we saw before 
this is essentially a normalization condition on the eigenstates. If X = {s}, then the 2|X|e(X) 
and 2e(X) terms in /(X) cancel. For such an X, e(X — m) is zero except when m = s, in 
which case it is 1. So for such X, /(X) = becomes 

e. = 2 h{Y)e{Z). (78) 

y,z,ieZ:(y+i)Az={s} 



For X with |X| 1 (including the case of X = 0), the equation /(X) = can be written as 

1 



e(X) 



2(|X|-1) 



-2 M^)e(^)+Eeme(X-m) 

Y,Z,j^Z:{Y+j)/\Z=X rn 



(79) 
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We will show these equations have a solution by writing them as a fixed point equation. 
Consider the set of variables 

e := {e{X) : |X| ^ 1} U {e, : s G A}. (80) 
Equations (78) and (79) form a fixed point equation for e 

F(e) = e. (81) 

Let us introduce the norm 

l|e||:=EKI(H^r'^' + 2 E |e(X)|||X|-l|(|6|M)-oW. (82) 

s X:\X\^1 

The factor of 2 is included merely for convenience. Here wo{X) denotes the number of bonds 

in the smallest connected set of bonds containing X and the origin, and |s| is the norm of 
the site s. Note that \s\ = Wo{{s}). li X — XiA ■ ■ ■ AX„, then 

WoiX) < WoiXi) + ■■■ + WoiXr,). (83) 

We prove that the fixed point equation for e has a solution by using the contraction mapping 
theorem as we did in the previous section. We must show that there is a 5 > such that 

||F(e) -F(e')|| < -||e-e'|| for | |e| |, | |e'| | < 5, (84) 

\\F{e)\\<6 for ||e|| < 5. (85) 
To verify (84), we use (78) and (79) to see that 

lira-F(e')|| 



< E (|e|^)""'°^^^ 

X:\X\j^l 



2 E \h{YMZ) - e'{Z)] I + E MX - s) - e',e\X - 

Y,Z,jGZ:{Y+j)AZ=X s 

+2E(k|M)-l'l E \h{YMZ)-e'{Z)]\. 

I Y,Z,jeZ:{Y+j)AZ={l} 



For a site /, we have Wo{{l}) = so the above 

= 2 E (|e|M)-"'°«^+^')^^)|/i(y)[e(Z) -e'(Z)]| 

Y,z,jez 

+ E E (|6|M)-oW|e,e(X-s)-e;e'(X-5)|. (87) 

s X:\X\^1 

We claim that for j G Z we have wq{{Y + j)AZ) < wq(Y) + wq{Z). To see this, let A and B be 
connected sets of bonds containing the origin such that Y <Z A, Z <Z B and 1^41 = wo{Y) and 



18 



\B\ — wo{Z). Since the origin is in ^4, the site j isin A + j. The site j is also in Z and so is in 

B. Since each of A -\- j and B is connected and they both contain j, their union is connected. 
So if we define C = {A + j) U B, then C is a connected set of bonds which contains the origin 
and {Y + j)AZ. So 

wo{{Y + j)AZ) < \C\ <\A + j\ + \B\^ \A\ + \B\ = wo{Y) + wo{Z). (88) 

We also have Wo{X) < Wo{X — s) + \s\. The sum over j in (87) then produces a factor of \Z\. 
So (87) is 

< 2^(|e|M)-^°('^)-"'°(^)|Z||/i(F)[e(Z) - e'{Z)]\ 

Y,Z 

+ E E (|6|M)--(^-)-N|e,e(X-5)-ey(X-5)|. (89) 

s X:\X\f^l 

Since e{Z) = e'{Z) if \Z\ = 1, and \Z\ <2\\Z\-1\ if \Z\ ^ 1, we have 
^(|e|M)-"'°(^)|Z||e(Z)-e'(Z)| < ^ (|e|M)-"'°(^)2||Z| - l||e(Z) - e'(Z)| < ||e - e'||. (90) 

So the first of the two terms in (89) is bounded by 

2$:(|6|M)-o(^)|MF)|||e-e'||. (91) 

Y 

After the change of variables X — > X + s, the second term in (89) equals 

E E (|e|M)--W-W|e,e(X)-e;e'(X)| 

<E E (|6|M)-°W-H[|e,e(X)-e,e'(X)| + |e,e'(X)-ey(X)|] 

s X:\X\i^\ 

< l|e-e'|| + ^||e-e'|| ||e'|| < (5||e-e'||, (92) 

where we have used the triangle inequality and the fact that ||e|| and ||e'|| are no greater than 
b. Combining (91) and (92), we have 

||F(e) -F(e')|| < i^||e-e'||, (93) 

where 

X = 2^(|e|M)-"'°(^)|/i(y)| + 5. (94) 

Y 

From (66) and (83) 

OO 1 

^(|e|M)-°(^)|/i(y)| < ^ - (|e|M)-o(^^)-— °(^")|^(yO • ■ ■ g{Y^)\ = ell^H - 1. 

y n=i ^- yi30,...,y„90 

(95) 
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So 

X < 2(ell»ll - 1) + 5. (96) 

If S and e are small enough, then K <l/2. 

To prove (85), we use (78) and (79) to compute F{0). Note that e = means that e{X) — 
for all X except X = {0}. We always have e({0}) = 1. Letting e' denote F{0), we have 

eH2 E h{Y) (97) 

Y:YA{0}={1} 

and for X with |X| 7^ 1 (including the case of X = 0), 

I -"-J Y:Y A{0}=X 

Thus 

||i^(0)|| < 2^(|e|M)-^°('*')|/i(y)| < 2(ell^ll - 1). (99) 

Y 

Hence, for ||e|| < 6 and e small enough that 2(e"^" — 1) < 6/2, 

\\F(e)\\ < \\F(e) - F(0)|| + ||F(0)|| < 1| |e| | + | |F(0)| | < 5. (100) 

This completes the proof that the fixed point equation has a solution. 

We still must prove that for each k, the eigenstate we have constructed has the lowest 
eigenvalue in the subspacc of momentum k if k ^ and has the next to lowest eigenvalue if 
k = 0. When e = we know that this is true. So if we can prove that within the subspace of 
momentum k, our eigenstate is never degenerate for |e| < eo, then we are done. We will use 
the local uniqueness of solutions to the fixed point equation to prove this. 

Let '0fc(cr) continue to denote the eigenfunctions constructed above with eigenvalues Ei{k). 
Suppose that for some ko there is another eigenfunction ijj'{a) with momentum ko and eigenvalue 
Ei{ko). We can write it in the form 

i^'(a)^^(a) Y: c'{X)<I>x,M- (101) 
rk{a) = M^) (102) 

^'^^""^^ l + ac'({0}) ■ ^ ^ 

(We need only consider small a, so we are not dividing by zero in the above.) Then for all 
small CK, ipki^) is eigenstate of H with eigenvalue Ei{k) and momentum k for all k. These 
states can be written in the form 

^W^) = ^{a) Yl Ca{X, k) cl>xM, (104) 



For small a we define 
for k ^ ko and 
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where Ca{X, k) is just c{X, k) for k ^ ko and 



Ca{X, ko) 



c{X, kp) + ac\X, kg) 
l + ad{{0}) 



(105) 



As before we can write this in the form 



(106) 



X 



where 



Ca{X, k) 



(107) 



n 



and 



ea{X) 



E 



e«(Y,n). 



(108) 



YeXM:Y+n=X 



We have normahzed things so that ea({0}) = 1 and ea({/}) = for / 7^ 0. It follows from what 
we have done before that for small a, e^^X) is a solution of the fixed point equation (75). As 
q; ^ 0, ea{X) — i> e(X). However, our solution of the fixed point equation is locally unique 
by the contraction mapping theorem. Thus ea{X) — e{X) for small a. This implies that the 
second cigenfunction ■0' is in fact just a multiple of ipko- This completes the proof that the 
eigenfunctions we have constructed are indeed the lowest excited states. 

The existence of the infinite volume limit and the convergence of to follow from the 
estimates of this section and the arguments of the previous section. In the fixed point equation, 
(78) and (79), h{Y) depends on the volume. However, the difference between it and its infinite 
volume limit is exponentially small in L. The extension of e{X) to a function e'{X) that is 
defined for all X is simpler than the definition of g'{X) in the previous section. We can simply 
define e'{X) to be e{X) if X is in A and otherwise. The bounds on the Cg follow from our 
choice of norm (82) which implies that 



As an application of the theorem, we show that the minimum of the dispersion relation 
occurs at A; = (tt, • • • , tt) if e > and at A; = (0, • • • , 0) if e < 0. We can use our expansion for 
the dispersion relation to compute the order e contribution to the dispersion relation. Consider 
eq. (78) for the coefficients e^. Recall that e({0}) = 1. Otherwise, e{Z) is always at least order 
e. For all Y, h{Y) is at least order e. So can get an order e contribution from the sum in (78) 
only for Z — {0}. In this case j must be zero. The constraint on the sum is yA{0} = {s}, 
which is equivalent to y = {0}A{s}. Prom (21) we see that the only Vs for which g{Y) has 
a nonzero contribution at first order in e are the sets containing two sites which are nearest 
neighbors. For such sets g{Y) = e/2 + O(e^). Thus from (66) we see that h{Y) has a nonzero 
contribution at first order in e only if F is a set consisting of two nearest neighbors, one of 




(109) 



and so 



65 is of order |e|l*l. This completes the proof of Theorem 2. | 
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which is the origin. For such sets h{Y) — e/2 + O(e^). Since Y — {0}A{s}, has a first 
order contribution only if s is a nearest neighbor of 0, i.e., |s| = 1. In this case eq. (78) gives 
e, = e + 0(e2). Thus 

d 

A{k) = 2 + e ^ e'''' + ^ e'.e'''' = 2 + 2e ^ cos{ki) + ^ e'.e''''. (110) 

s:\s\=l s&Zd i=l seZ'' 

where the are all of order e^, and ki are the components of the vector k. 

We consider the case of e > 0. The argument in the other case is analogous. Obviously, 
2 + 2e J2i=i cos{ki) has its minimum at k = {%, ■ ■ ■ ,%). To show this is true for A(/c) we consider 
two cases. It is easy to see that /c = (tt, • • • , tt) is a critical point of A(A;) and by computing the 
Hessian that this critical point must be a local minimum. So there is a 5 > so that A(/c) is 
greater than its value at (tt, ■ ■ ■ , vr) provided \ki — it\ < 6 for all i. If any ki has |A;j — 7r| > S, 
then 2 + 2e X^^Li cos(A;j) is greater than its minimum value by at least 0(5^). The difference 
between 2 + 2€Y,i=i cos{ki) and A(A;) is O(e^), so \ki — tt] > 6 implies A{k) is greater than its 
value at A; = (tt, • • • , tt) provided e is small compared to 5. 
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